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Abstract. The Bruggeman formalism is implemented to estimate the refractive index of an 
isotropic, dielectric, homogenized composite medium (HCM). Invoking the well-known Hashin- 
Shtrikman bounds, we demonstrate that the group velocity in certain HCMs can exceed the group 
velocities in their component materials. Such HCMs should therefore be considered as metamate- 
rials. 

PACS numbers: 41.20.Jb, 42.25.Dd, 83.80.Ab 

1 Introduction 

By definition, metamaterials exhibit behavior which (i) either their component materials do not 
exhibit (ii) or is enhanced relative to exhibition in the component materials Many types of 
metamaterials may be conceptualized through the process of homogenization 121 011], paving the 
way for their realization. For example, a homogenized composite medium (HCM) may be envisaged 
which supports the propagation of a Voigt wave (which is a planar wave whose amplitude varies 
linearly with propagation distance), although such waves cannot propagate through its component 
materials jSlini- 

In this communication, we explore the enhancement of group velocity which may be achieved 
through homogenization. S0lna and Milton recently considered this issue, by estimating the relative 
permittivity of a HCM as the volume-weighted sum of the relative permittivities of the component 
materials 0. But that estimation is applicable only for planar composite materials such as super- 
lattices of thin films, and not to the more commonly encountered particulate composite materials 
[21IS1I5- III the following analysis, we implement the well-established Bruggeman formalism [1] 
to calculate the effective refractive index of an isotropic dielectric HCM. Thereby, we demonstrate 
that metamaterials which support group velocities exceeding those in their component materials 
may be realized as particulate composite materials. 
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2 Analysis 



Consider a composite material containing materials labeled a and b, with refractive indexes 
and Hf), respectively. The component materials are envisioned as random distributions of spherical 
particles. Provided that the diameters of these particles are small compared with electromagnetic 
wavelengths, homogenization techniques may be applied to estimate the effective refractive index 
of the HCM. 

In particular, the well-established Bruggeman homogenization formalism [21 EI — which may 
be rigorously derived from the strong-permittivity-fluctuation theory TKSl^ M03 — leads to the 
equation 
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whose solution yields nsr as the estimated refractive index of the HCM. Here, fa and fb = ^ — fa 
are the volume fractions of the component materials. In the following, both component materials 
are assumed to have negligible dissipation in the frequency range of interest. 



The group velocity of a wavepacket propagating through the HCM is given as |lUj 



VBr = 



(2) 



where vbt is evaluated at the angular frequency uj = u){kavg), with kavg being the average wavenum- 
ber of the wavepacket, and c is the speed of light in free space. Similarly, the respective group 
velocities in component materials a and b are given by 



(i =a,b). 
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We proceed to establish upper and lower bounds on VBr, in terms of Ua and Uf,- In particular, we 
demonstrate that the inequalities 



VBr > Ve , 



{£ = a, b) 



can be satisfied for certain values of > 1, rif, > 1, > 0, and ^ > 0. 
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Differentiation of both sides of Eq. with respect to uj yields 
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Upper and lower bounds on 6a and 5b may be established by exploiting the Hashin-Shtrikman 
bounds riL and nu on nsr i.e., 

ni < riBr < nu, 



where 
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Combining Eqs. (jEI)^©, we get 

Pi < 5i < Ke, 

where 
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Thus, we have 
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and the group velocity in the HCM is accordingly bounded as 



VL < VBr < VU , 



with 
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If the inequalities 
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hold for certain component materials, then the inequalities (0} are automatically satisfied. 
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The inequalities ()15() reduce to the particularly simple inequality 

nu + ^-r^ {Ka + Kt- 1) + Kb {Ua - rife) < Ua, (16) 

if Va = Vh- The conditions 

nu + Kb (Ua - Ub) < Ua 
Ka + Kb-l > 

are satisfied, for example, by Ua = 3, n;, = 1.2 and fa = 0.9. Thus, the inequality (|16j) holds, 
provided that the dispersive term is sufficiently small. 




3 Numerical results 

Let us illustrate the phenomenon represented by the inequalities by means of a specific numerical 
example. Consider a particulate composite material at a particular value loq oi uj. At the chosen 



angular frequency, let = 5, n;, = 1.2, 



= 0.5M, and ^ 

LU=UJQ 



= 5.5/a;o. Significantly, 

material a has a high refractive index but low dispersion in the neighbourhood of Wq, whereas high 
dispersion in material b is combined with a low refractive index. 



The Bruggeman estimate of the refractive index of the HCM, namely UBr-, is plotted as a function 
of the volume fraction fa in figure 1. Also shown are the upper and lower Hashin-Shtrikman bounds, 
njj and ni, on ubt-, as well as the parameters 5a and 5b- The Bruggeman estimate adheres closely 
to the lower bound ni, at low values of fa, whereas at high values of fa the difference between ubt 
and its upper bound njj becomes marginal. The observed agreement between n^r and ni aX, low 
fa reflects the fact that the lower Hashin-Shtrikman bound is equivalent to the Maxwell Garnett 
estimate of the refractive index of the HCM arising from spherical particles of material a embedded 
in the host material h [3]. The Maxwell Garnett estimate is only valid then at low values of fa- As 
the volume fraction becomes increasingly small, the Bruggeman estimate {nsr) and the Maxwell 
Garnett estimate (low fa value of riL) converge on Ub- In a similar manner, the agreement between 
nsr and njj at high values of fa is indicative of the fact that the upper Hashin-Shtrikman bound 
is equivalent to the Maxwell Garnett estimate of the refractive index of the HCM arising from 
spherical particles made of material b embedded in host material a; the Maxwell Garnett estimate 
then holds only at high values of fa- In the limit fa 0, the coefficients 5a ^ and (5^ ^ 1; while 
5a ^ 1 and 5b ^ as fa ^ 1. 

The corresponding group velocities Va, Vb and VBr are plotted as functions of fa in figure 2. The 
upper and lower bounds on VBr as given by vu and vl, respectively, are also displayed. Clearly, we 
have VBr > Va and vbt > Vb for > 0.67. 

The inequalities @ hold only over a relatively small range of parameter values. For example, 
the phase space in which the inequalities (jU are satisfied is illustrated in figure 3 for = 5, 

0.5/tJo and fa = 0.8. With these relationships fixed for the component material a, we 
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that VBr > Va and vbt > Vb for 



(i) 1.17 <nb< 1.23 with ^ = S.Sl/cJo; and 
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< 5.57/a;o with = 1.2. 

4 Concluding remarks 

We conclude that the group velocity in an isotropic, dielectric, particulate composite material — as 
estimated via the Bruggeman homogenization formalism — can exceed the group velocities in its 
component materials. This metamaterial characteristic may be achieved through homogenizing (i) a 
component material a with high refractive index and low dispersion with (ii) a component material 
h with low refractive index and high dispersion. Neither anomalous dispersion nor an explicit 
frequency-dependent model of the refractive index (unlike Ref . [Jj ) is required to demonstrate this 
characteristic. 

Improved estimates of HCM group velocity may be achieved through the implementation of 
homogenization approaches which take into better account the distributional statististics of the 
component materials (e.g., the strong-permittivity-fluctuation theory approach TK81j M03). In 
particular, the effects of coherent scattering losses — which are neglected in the present study — 
may well result in a moderation of the group velocity. Such studies are currently being undertaken, 
especially in light of the recent emergence of metamaterials wherein the phase velocity and the 
time-averaged Poynting vector are oppositely directed LM^03. 
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List of Figure Captions 



Fig. 1. The estimated refractive index nsr (solid line), the upper and lower Hashin-Shtrikman 
bounds on n^r (broken dashed lines, labeled as nu and n^), and the coefficients 6a and 6b (dashed 
lines), all plotted as functions of the volume fraction /„, when Ua = 5 and = 1.2. 

Fig. 2. The estimated group velocity vbt (solid line) and its upper and lower bounds (broken 
dashed lines, labeled as vu and vl), along with the group velocities Va and vi, (broken dashed lines) 
in the component materials, plotted as functions of the volume fraction /„, when Ua = 5, = 1.2, 
0.5/^0 and ^ = 5.5/a;o. All group velocities are normalized with respect to c. 



du! 



u>=u>o 
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Fig. 3. The shaded region indicates the portion of the a-/3 phase space wherein vbt > Va and 
VBr > Vb'i here, a = ^ and P = (^/^)\ ■ This region was demarcated for Ua = 5, 



dw 



= 0.5/^0 and = 0.8. 

U1=U!0 
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Figure 1: The estimated refractive index UBr (solid line), the upper and lower Hashin-Shtrikman 
bounds on nsr (broken dashed lines, labeled as nu and n^), and the coefficients 5a and 6h (dashed 
lines), all plotted as functions of the volume fraction /„ when ria = 5 and = 1.2. 
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Figure 2: The estimated group velocity vbt (solid line) and its upper and lower bounds (broken 
dashed lines, labeled as vij and u^), along with the group velocities Va and (broken dashed lines) 
in the component materials, plotted as functions of the volume fraction /„, when = 5, nj, = 1.2, 

S.S/wq. All group velocities are normalized with respect to c. 
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Figure 3: The shaded region indicates the portion of the a- [3 phase space wherein v^r > Va and 
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